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ABSTRACT 

The central star of the planetary nebula PRTM 1 (PN G243. 8-37.1) was previously found to be variable by M. Pena and colleagues. 
As part of a larger programme aimed towards finding post common-envelope binary central stars we have monitored the central star 
of PRTM 1 spectroscopic ally and photometrically for signs of variability. Over a period of ~3 months we find minimal radial velocity 
(<10 km s _1 ) and photometric (< 0.2 mag) variability. The data suggest a close binary nucleus can be ruled out at all but the lowest 
orbital inclinations, especially considering the spherical morphology of the nebula which we reveal for the first time. Although the 
current data strongly support the single star hypothesis, the true nature of the central star of PRTM 1 remains enigmatic and will 
require further radial velocity monitoring at higher resolution to rule out a close binary. If in the odd case that it is a close binary, it 
would be the first such case in a spherical planetary nebula, in contradiction to current thinking. 
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1. Introduction 

At least 17 ± 5% of planetary nebulae (PNe) have post common- 
envelope (CE) central stars with orbital periods less than ~1 
day ([Miszalski et al. , 2009a ). After a long and protracted debate 
(Balick & Frank 2002), binary central stars are now emerging as 
the preferred mechanism to shape PNe (Soker 2006; Nordhaus 
& Blackman 2006; De Marco 2009 and references therein), long 
after their predicted existence and subsequent discovery in the 
mid-70s (Paczynski 1976; Bond, Liller & Mannery 1978). In or- 
der to understand the full influence of binaries on the formation 
and evolution of PNe a large and representative sample of bi- 
naries must be identified. Since the discovery of UU Sge as an 
eclipsing nucleus of the planetary nebula Abell 63 (Bond 1976, 
Bond et al. 1978), there are at least 40 close binary central stars 
now known, the overwhelming majority of which were found 
via photometric monitoring to reveal periodic variations due to 
irradiation, ellipsoidal modulation and eclipses (e.g. Miszalski 
et al. 2008, 2009a, 2011a, 2011b, 2011c; Corradi et al. 2011; 
Santander-Garcia et al. 2012, in prep.). Many of these show a 
strong irradiation effect, due to the atmosphere of a K- or M- 
type dwarf companion being heated by the white dwarf primary, 
an effect that becomes undetectable once the orbital period ex- 
ceeds a few days (De Marco, Hill wig & Smith 2008). A strong 
selection effect therefore biases our current knowledge of binary 
central stars which may reside within a much larger fraction of 
the PN population (Moe & De Marco 2006, 2012). 

Radial velocity (RV) monitoring of central stars is an under- 
utilised methodology which has the capability to detect binaries 
that lack a favourable combination of orbital period and incli- 
nation to show detectable photometric variations. Boffin et al. 



* Based on observations made with the Very Large Telescope at 
Paranal Observatory under programmes 088.D-0750(B) and 078.D- 
0033(A). 



(2012) powerfully demonstrated the potential of RV monitor- 
ing by discovering the double-degenerate nucleus of Fleming 1 
which showed a flat lightcurve and no irradiated emission lines 
(e.g. Miszalski et al. 2011b, 2011c). Previous RV monitoring 
surveys have unfortunately not been as successful with no rou- 
tine discoveries made so far. Mendez (1989) estimated 15% of 
his monitored sample to show RV changes but none were pe- 
riodic. High levels of variability were also found by De Marco 
et al. (2004), Sorensen & Pollacco (2004) and Afsar & Bond 
(2005), but again no strictly periodic binaries were found (see 
also Frew et al. 2010). Although this could in principle partly be 
due to intrinsic, i.e. wind, variability, this can perhaps be mostly 
explained by the sparse sampling employed in these surveys, 
whereby spectra were often taken days or months apart, severely 
reducing their sensitivity to short periods. Improved sampling as 
in the Boffin et al. (2012) study is paramount to improving the 
utility of RV monitoring of central stars. 

The selection of Fleming 1 for RV monitoring by Boffin et al. 
(2012) was based on the extraordinary point- symmetric jets of 
Fleming 1 (e.g. Lopez et al. 1993) and the unusual spectroscopic 
type of its central star. Its spectroscopic appearance matched that 
of the variable central star of the PN PRTM 1 (PN G243. 8-37.1, 
Pena et al. 1990 and Pena & Ruiz 1998, hereafter PRTM90 
and PR98, respectively). If this unusual spectral type were to 
be a signature of a close binary nucleus, then PRTM 1 may also 
be expected to have one. Indeed, this possibility was raised by 
PR98 to explain the spectroscopic and photometric variability of 
PRTM 1 . Following the successful monitoring of Fleming 1 , we 
decided to test this hypothesis with time- series spectroscopy of 
the central star of PRTM 1 over 3 months. Here we report on 
the results of our campaign and present the first morphological 
study of the surrounding nebula. 

This paper is structured as follows. Section [2] presents a 
very brief summary of previous work on this remarkable ob- 
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Fig. 1. Rectified average FORS2 spectrum of the central star of PRTM1. The nebular contribution has been mostly removed. 



ject. Section |3]presents our spectroscopic and photometric study 
of the central star using the FORS2 instrument at the VLT. 
Section [4] describes our imaging and kinematics of the nebula. 
A discussion of our results is in Sect.[5]and Sect. [6] concludes. 



2. Previous work on PRTM 1 

Both the central star and nebula of PRTM 1 have been stud- 
ied previously by numerous authors. PRTM 1 was discovered 
by PRTM90 from objective prism survey plates and resembles 
other highly ionised Halo PNe (e.g. NGC 4361). The nebula 
spectrum is highly ionised and demonstrates only a slightly sub- 
solar abundance pattern, unusual for a Halo PN at a distance 
of 5 ± 2 kpc (PRTM90). PR98 reported the central star to have 
an O(H) spectral type with T e s = 80 kK and log g=5.2 based on 
analysis of the central star by Mendez (1991)QPR98 found vari- 
ability in the central star based on low-resolution spectroscopy 
which shows a blue continuum, deep H I and He II absorption 
lines, as well as faint and narrow C IV AA5801, 5812 emission 
lines. The lack of C III /15696 emission rules out a [WC] classi- 
fication, while non-detection of O VI /L1381 1-34 emission rules 
out a [WO] classification (Crowther, De Marco & Barlow 1998). 

Feibelman (1998) found that IUE observations obtained in 
1988 and in 1990 showed significant differences, the stellar flux 
increasing by a factor of ~3, while the He II Al 640 emission 
line was 1.5 times higher. PR98 also noticed important changes 
in the spectra of the central star, both in the continuum and ab- 
sorption line profiles, when comparing spectra obtained on 28 
October 1990, 1 February 1998 and 22 March 1998. Moreover, 
PR98 noted large variations in the stellar flux of the PN, with 
three deep minima that occurred on 20-22 November 1988, 10 
August 1990, and 28 October 1990. The largest difference in 
brightness corresponded to AV~1 mag. PR98 concluded that the 



1 Mendez (1991) referenced Mendez et al. (1988) for this analysis. 
However, we did not find PRTM 1 to be included in this work. 



stellar continuum changed by a factor 2-3 on short timescales 
(weeks) and 20%-30% variations occurring within days. 

PR98 proposed 3 possible scenarios to explain these varia- 
tions: (a) a close binary system with a period of a few hours to 
a few days that shows eclipses or an irradiation effect, (b) non- 
radial pulsations in a very hot high-gravity pre- white dwarf, or 
(c) stellar wind variations. Our observational campaign aimed to 
test their preferred scenario of a close binary central star (a). 



3. The central star 

PRTM 1 was observed with the FORS2 instrument (Appenzeller 
et al. 1998) under the Very Large Telescope (VLT) service mode 
programme 088.D-0750(B). Table [T] lists the 15 useful spectra 
taken over the course of 84 days from 11 October 2011 to 3 
January 2012. A 0.5" slit and the 1400V volume-phase holo- 
graphic grating were used to give a wavelength range of 4597- 
5901 A with a dispersion of 0.32 A pixel -1 and a spectral reso- 
lution measured from the arc spectrum of 1 .22 A (full width at 
half-maximum, FWHM). Exposure times mostly varied between 
900-1500 s resulting in an average signal-to-noise ratio (S/N) of 
105 near 5100 A. Basic reductions of the data were performed 
using the ESO FORS pipeline and one-dimensional spectra were 
optimally extracted using midas. The mean spectrum of the cen- 
tral star of PRTM1 is shown in Fig.[T] Since it is not possible to 
take contemporaneous arc exposures with FORS 2, the error bar 
on the heliocentric radial velocity is rather high, and we there- 
fore do not attempt at measuring it (but see Sect. 4.2). 

In view of the results of PR98 who found large variations in 
the spectroscopic behaviour of the central star, we compared our 
higher resolution spectra against theirs. No major changes were 
seen over the full 84 day monitoring period. We also measured 
the equivalent width of the various absorption lines and did not 
see any significant change. The absorption line at 5114.06 A, 
most likely due to O V, has a rather constant equivalent width of 
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Table 1. Spectroscopic observations and measured radial veloci- 
ties of the central star of PRTM 1. The table shows, respectively, 
the modified julian date (MJD) at the midpoint of each observa- 
tion, the exposure time, the attained signal-to-noise ratio and the 
measured radial velocity (see text). We estimate a typical error 
of 3 km s" 1 on the latter. 
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Fig. 2. Photometric and spectroscopic measurements of PRTM 1 
during our 84-day campaign. The top panel shows the estimated 
Rf, Bj, and In magnitudes as a function of time, while the bot- 
tom panel shows the radial velocity of the O V 5114.06 absorp- 
tion feature in the rest- frame of the nebula (see text). 



85 mA, the l<x variation being of 7.7 mA, i.e. a change of ~9%, 
well below the factor 2-3 found by PR98. 

The FORS2 spectra were converted to a logarithmic wave- 
length scale before shifting each spectrum to match the systemic 
nebula RV based on the strong [O III] emission lines. Gaussian 
profile fitting of the nearby O V A5 114.06 absorption line was 
used to measure the central star's RV with respect to the nebula 
RV. Table [T] lists the derived RVs together with the S/N reached. 
Similar RV variations were found from other stellar features (e.g. 
He II ,15412, C IV A5801, ^5812) which all moved in phase with 
O V. From the signal-to-noise of the spectra and the spectral res- 



olution, we estimate, from an extensive set of Monte Carlo simu- 
lations, an error of 3 km/s per data point - or slightly smaller for 
the best of our spectra. This error encompasses both the (small) 
error related to the position of the nebula emission line and that 
related to the position of the O V line. The slit used was 0.5 
"wide and was therefore always smaller than the seeing. 

Figure [2] shows the RV measurements over time. The RV is 
apparently changing with time, spread over a range from -14.9 
to +6.0 km s" 1 , with variations occurring on a time scale of a 
few days. The mean value is -3.7 km s" 1 and we measure an 
intrinsic scatter of 6.6 km s" 1 , much larger than the individual 
errors. The figure seems to indicate a slow upward trend, but this 
may simply be due to the limited sampling. To quantify whether 
this trend is real, we have compared the residuals when removing 
or not the linear trend identified in Figure[2j the^ 2 changes from 
47 to 29, and a F-test applied to these number reveals that this is 
not significant (probability = 0.38). We therefore conclude that 
this trend is most likely spurious. 

During the campaign a variety of acquisition images were 
also taken to monitor any photometric variability and to study 
the morphology of the nebula (Sect. [4]). Table|2]gives the full list 
of acquisition images taken. The observations were performed 
using the filters 011+44 (central wavelength, A c = 371.7 nm, 
width, W = 13 nm), OIII+50 (A c = 500.1 nm, W = 5.7 nm), 
H_Alpha+83 (A c = 656.3 nm, W Q = 6.1 nm) and LBESS+77 
(A c = 768.0 nm, W = 138.0 nm). An approximate absolute 
scale to our photometry was determined by using the catalogue 
photometry of the SuperCOSMOS Sky Survey (Hambly et al. 
2001), namely I N for LBESS+77, R F for H_Alpha+83 and Bj 
for OIII+50. Differential photometry was performed using more 
than 10 field stars in all filters, using the same methodology as 
that used by Miszalski et al. (201 lc), i.e. using sextractor and 
an aperture of radius of 3" (-1.5 the worse seeing). We also 
made an analysis using standard aperture photometry and ob- 
tained similar results. There were insufficient stars in the [O II] 
frames to derive an appropriate calibration. Note that because of 
the large aperture used, our values partly include a nebula con- 
tribution (affecting the Bj and R F magnitudes the most). These 
magnitudes are included in Tableland displayed alongside the 
RVs in Fig. [2] There is no evidence for a change in central star 
brightness exceeding ~0.2 mag during the time of our spectro- 
scopic monitoring. 



4. The nebula 

4.1. Morphology 

Figure [3] shows the first high-resolution images taken of PRT M 1 
featuring the observations with the best seeing in Table [2p] The 
data show a filamentary inner nebula 14-1 5" across surrounded 
by a fainter shell measuring 20.5 x 21.3 arcsec. At the edges 
of the apparent major axis along PA~135° there seems to be 
two thin arcs in the [O III] image. The nebula is typical of 
multiple- shell PNe (Chu, Jacoby & Arendt 1987) and the some- 
what 'dented' appearance of the inner nebula is unusual but also 
seen in other PNe (e.g. NGC 2867, Schwarz, Corradi & Melnick 
1992). Given the very high degree of nebula ionisation there 
are no low-ionisation structures present and no faint features 
located outside the fainter shell. Together all these properties 
are inconsistent with the nebula morphologies of post common- 
envelope PNe (Miszalski et al. 2009b). It is interesting to note 



2 Specifically data taken on the nights of MJD=55882 ([O II]), 
MJD=55905 ([O III]) and MJD=55906, 55920 (Ha). 
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Table 2. Log of acquisition images and derived magnitudes for PRTM 1 . The magnitudes include nebula and central star in the 
aperture. 
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Fig.3. VLT FORS2 images of PRTM1. (a) Ha, (b) [O III], (c) [O III] with the ARGUS 11.5 x 7.3 arcsec 2 footprint (see Sect.|4.2|), 
(d) [O II], (e) Ha divided by [O III] (Ha excess appears light and [O III] dark), and (f) Colour-composite made from Ha (red), 
[O III] (green) and [O II] (blue). The images measure 90 x 90 arcsec 2 with North up and East to left. 



that several of the nebulae surrounding GW Vir pulsating central 
stars strongly resemble PRTM 1 including NGC 7094, Abell 43, 
NGC 2867 and especially NGC 1501 (Fig.g]). 



4.2. Kinematics 

To study the kinematics of the nebula of PRTM 1, we retrieved 
high resolution observations from the ESO archive. The PN was 
observed on 15 January 2007 using the FLAMES facility of the 
VLT (Pasquini et al. 2002) under programme 078.D-0033(A) 
(PI: D. Schonberner, see Schonberner et al. 2010). FLAMES 
was operated in its ARGUS mode, namely an integral field unit 
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Fig. 4. NOAO Hot image of NGC 1501 (Credit: Jay Gallagher 
U. Wisconsin/WIYN/NOAO/NSF). The diameter of the nebula 
is -60" across and North is up with East to left. Note the strong 
resemblance to PRTM 1 (Fig. [3). 

(IFU) of 22 x 14 microlenses, each 0.52" a side, giving a foot- 
print of 11.5 x 7.3 arcsec 2 (see Fig. [3^). The HR8 grating was 
selected for the GIRAFFE spectrograph providing a resolving 
power R = 32, 000 around the [O III] A5001 emission line. Two 
1800-s exposures were retrieved from the ESO archive and re- 
duced using the ESO pipeline before being averaged using stan- 
dard starlink routines. The seeing during the observations var- 
ied between 1.0-1.5". 

Figure [5] shows Gaussian fits to an extracted spectrum equiv- 
alent to a 2 spaxel-wide (1.04") longslit spectrum through the 
CSPN at a position angle of 90°. The velocity profile shows lit- 
tle deviation from that expected from a spherical nebula, strongly 
indicating that PRTM 1 is roughly spherical with an expansion 
velocity of 30.0 ±5.0 km s" 1 and a heliocentric systemic ra- 
dial velocity of 12.5 ±5.0 km s" 1 . The relatively large error 
bars reflect turbulence in the nebula. Assuming a distance of 
5 ± 2 kpc (PRTM90), the kinematic age of the -0.5 pc diam- 
eter nebula corresponds to 5500^^ years including both errors 
in the expansion velocity and distance. Both the kinematics and 
morphology of PRTM 1 are consistent with those of NGC 1501 
(Sabbadin & Hamzaoglu 1982). 

5. Discussion 

In combination with the observations of Pena and colleagues, 
the observations presented in this paper allow us to refine the 
possible nature of the central star. We discuss three possible sce- 
narios: a post common-envelope binary, a wider binary with an 
orbital period of weeks or months to years and a single central 
star. 

5.1. A close binary central star? 

A thorough periodogram search was performed, both on the ob- 
served RV time- series and the corrected series where we re- 
moved the spurious long-term trend. Although some periods 
phased up reasonably well (e.g. P~3 hours), these always had at 
least 2 uncooperative RVs that could not be explained away - i.e. 
they were more than 5-cr outliers. In any case, these periods had 
too low a S/N in the periodogram, meaning that no clean periods 
were found. If the RV variations we have observed were due to 
the presence of a companion, the orbital period would have to be 




-40 -20 20 40 60 

Heliocentric radial velocity (km s~ 1 ) 



Fig. 5. VLT ARGUS velocity profile of PRTM 1 from a 1.04" 
longslit trace, overlaid with the profile expected for a spherical 
nebula at a heliocentric systemic velocity of 12.5 km s" 1 . Note 
that for clarity and ease of comparison to a spherical model only 
one PA is presented, but the kinematics were examined across 
the entire nebula and are entirely consistent with the PA pre- 
sented. 



smaller than around a few days, to account for the wide changes 
seen in this time range. The semi-amplitude we measure would 
be on the order of 10 km s" 1 , which given reasonable values 
for the masses of the components would mean that the system 
is seen almost face-on. We can in fact try to be more quantita- 
tive. Using the definition of the RV semi- amplitude in a circular 
Keplerian orbit, we can define a relation between the orbital pe- 
riod and the inclination of the system in the plane of the sky, 
assuming the masses of the components. This is shown in Fig.[6j 
where we have taken 0.6 M Q as the mass of the white dwarf, 
and considered companions of masses 0.3, 0.9 and 1.5 M , re- 
spectively. The last value is probably unrealistically too high as 
PRTM1 is a Halo PN, whose progenitor must have had an initial 
mass well below this value and there is no reason to think that 
the secondary became more massive than the primary. Not sur- 
prisingly given the very low semi-amplitude measured, for most 
of the values of the inclination angles, the orbital period would 
have to be rather large, above 100 d (for M2=0.3 M©) or even 
1000 d (for M2=1.5 M©). As mentioned above, the variations 
we detect are on a much smaller time scale than those detected 
by PR98 and Feibelman (1998) and therefore could not be due to 
such a long period. Only if the inclination angle were very small 
would the orbital period be below a few days for the whole range 
of masses. The almost perfectly round shape of the nebula un- 
fortunately does not provide any clue on the possible presence 
of a torus or ring, and thus no clue to obtain the inclination of 
the orbital plane if it exists (see e.g. Jones et al. 2012; Tyndall et 
al. 2012). 

None of the spectra obtained here or in PR98 show emis- 
sion lines of C III and N III that would come from an irradi- 
ated atmosphere of a main-sequence companion (e.g. Miszalski 
et al. 2011c). Nor do we find any spectral features that can be 
attributed to a companion, all stellar features move in lockstep 
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Fig. 6. Possible values of the orbital period as a function of the 
inclination angle for 3 values of the secondary mass. The inset 
shows a zoom for the smallest inclination angles and periods. 



with each other. The roughly spherical nebula morphology and 
kinematics (Sect. [4} are also incompatible with a close binary 
on theoretical grounds (e.g. Soker 1997) and the observed mor- 
phologies of nebulae around currently known close binaries (e.g. 
Miszalski et al. 2009b; Miszalski 2012). While we cannot rule 
out a very low inclination close binary based on the relatively 
small number of observations currently in hand, it seems un- 
likely that future observations would find one. If there were a 
close companion it would probably be a low mass white dwarf. 

5.2. A wide binary central star? 

If the period were weeks to months we may have expected to see 
more coherent radial velocity variations over the ~3 months we 
observed the central star. Moreover, the relatively large level of 
spread in the RVs of cr ~10 km/s happening on a timescale of 
a few days appears to be real. If some larger periodicity would 
have to be invoked, it would have to be in addition to some other 
intrinsic variation. As we indicated above, the RV data shows 
some apparent upward trend. If we discard the possibility that 
this is simply due to small number statistics and our sampling 
- as we have shown above that it probably is - we may infer 
that the mean velocity increased by 15.7 km s" 1 over the 84-day 
period. Thus the orbital period in this case would be more than 
twice this, and the radial velocity amplitude should be of the or- 
der of 7-10 km s" 1 . The longer the period, the larger the radial 
velocity amplitude would have to be. As such amplitude is, how- 
ever, inversely proportional to a power of the period, too large a 
period would be impossible. With our values, we therefore es- 
timate that the largest possible period would be of -150-200 
days. 



5.3. A single central star? 

If PRTM 1 were to have a single central star, then we would have 
to explain the apparent RV variations as due to either pulsations 
or perhaps wind variability. We can probably rule out non-radial 
GW Vir type pulsations since the log g of 5.2 is ~ 2.0 dex too low 
for the r e ff = 80 kK temperature to place it in the expected in- 
stability domain (Fontaine & Brassard 2008). Furthermore, the 
O(H) spectral type of the central star of PRTM 1 is inconsis- 
tent with a pulsator since it cannot be classified as a PG1159 
(Napiwotzki & Schonberner 1995) or an early- [WC] central 
star (Crowther et al. 1998). Both of these spectral types domi- 
nate the sample of known GW Vir pulsating central stars (e.g. 
Ciardullo & Bond 1996). Wind variability would only be ex- 
pected if PRTM 1 were an Of(H) type central star (e.g. Mendez, 
Herrero & Manchado 1990). However, the spectrum clearly does 
not have this spectral type (Fig. [TJ. Another possibility is that 
secular evolutionary changes could explain the small amount of 
photometric variability constrained by our observations and the 
previously observed high levels of variability (see Sect. 1), as 
seen in some objects but on a much smaller scale (e.g. Bond 
2008). Such changes may be one possible explanation as to why 
we find no variability compared to previous authors. Finally, we 
re-emphasise the spherical morphology of the nebula that is con- 
sistent with a single central star. 

6. Conclusions 

We have presented in this paper a new study of the nebula and 
central star of PRTM 1 . Our imaging reveals a rather spherical 
morphology, and based on archival data, the kinematics are also 
consistent with a spherical nebula. We also conducted a spec- 
troscopic monitoring campaign of the central star of PRTM 1. 
Unlike previous authors, we do not find large variations in the 
continuum or spectral features over the full time range of our 
observations (84 days). We notice, however, some radial veloc- 
ity variations of the CSPN but were unable to find any coherent 
periodicity. At this stage, it is not possible to decide whether 
the CSPN is a binary or not, however the binary hypothesis 
seems very unlikely. If PRTM 1 were indeed a close binary, then 
it would be the first one discovered within a spherical nebula 
in contradiction with our understanding of how close binaries 
should shape PNe (e.g. Soker 1997). 
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